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ABSTRACT 

Context. It will soon become possible to directly link the most accurate radio reference frame with the Gaia optical reference frame 
using many common extragalactic objects. It is important to know the level of coincidence between the radio and optical positions of 
compact active galactic nuclei (AGNs). 

Aims. Using the best catalogues available at present, we investigate how many AGNs with significantly large optical-radio positional 
offsets exist as well as the possible causes of these offsets. 

Methods. We performed a case study by finding optical counterparts to the International Celestial Reference Frame (ICRF2) radio 
sources in the Sloan Digital Sky Survey (SDSS) Data Release 9 (DR9). The ICRF2 catalogue was used as a reference because 
the radio positions determined by Very Long Baseline Interferometry (VLB I) observations are about two orders of magnitude more 
accurate than the optical positions. 

Results. We find 1297 objects in common for ICRF2 and SDSS DR9. Statistical analysis of the optical-radio differences verifies that 
the SDSS DR9 positions are accurate to ~55 milliarcseconds (mas) in both right ascension and declination, with no systematic offset 
with respect to ICRF2. We find 5 1 sources (~4% of the sample) for which the positional offset exceeds 170 mas (~3cr). Astrophysical 
explanations must exist for the majority of these outliers. There are three known strong gravitational lenses among them. Dual AGNs 
or recoiling supermassive black holes may also be possible. 

Conclusions. The most accurate Gaia-VLSl reference frame link will require a careful selection of a common set of objects by 
eliminating the outliers. On the other hand, the significant optical-radio positional non-coincidences may offer a new tool for finding 
e.g. gravitational lenses or dual AGN candidates. Detailed follow-up radio interferometric and optical spectroscopic observations are 
encouraged to investigate the outlier sources found in this study. 
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1 . Introduction larly observed with Very Long Baseline Interferometry (VLB I) 

J-_ over a long period of time. 
>0 An ideal quasi-inertial reference system would be defined by 

point-like sources radiating in all wavebands of the electromag- Links between the ICRF2 and optical reference frames can 

O netic s P ectmm , providing a simple relationship between the po- be either direct or indirect. A direct connection would mean 

C<") sitions determined at different frequencies (Walter & Sovers that the positions of the radio AGNs are directly measured in 

2000). This would facilitate the cross-identification of celes- the optical band. At present, however, the most accurate optical 

tial objects in the radio, infrared, optical, ultraviolet, and high- astrometric catalogues cannot be linked directly to the ICRF2, 

. h energy bands. However, in practice, reference frames at differ- because their limiting magnitude significantly exceeds the de- 

^ ent frequencies are defined by different objects, since suitable tection limit of AGNs in the optical. For example, in the Sloan 

H sources usually radiate intensely only in particular wavebands Digital Sky Survey (SDSS) Quasar Catalogue (Paris et al. 2012), 

' (e.g. quasars can be bright in the radio but are relatively faint in m q >17 mag. On the other hand, the primary astrometric refer- 



the optical, while stellar objects in the Galaxy are usually bright ence used in the optical domain, the Hipparcos catalogue (Perry- 
in the optical but weak radio emitters at best). This makes it man et al. 1997) has a limiting magnitude of only V~ 12.4 mag. 
necessary to link the reference frames used in different wave- It is defined by bright stellar sources and could only be con- 
bands. The link between the optical and radio domains is partic- nected to the ICRF indirectly (Kovalevsky et al. 1997; Stone 
ularly important, since a dominant fraction of astronomical ob- 1998). Since no extragalactic radio source could be reliably de- 
servations is made in the optical, but the most accurate reference tected with Hipparcos, a variety of secondary methods had to be 
frame is currently realized in the radio. used, e.g. relative astrometric measurements of radio stars with 
The latest realization of the International Celestial Reference respect to nearby reference quasars with VLBI and connected- 
System (ICRS) at radio frequencies was adopted by the Interna- element radio interferometers, sensitive optical observations of 
tional Astronomical Union (IAU) in 2009. The second version quasars relative to Hipparcos stars, and a comparison of VLBI- 
of the International Celestial Reference Frame (ICRF2, Fey et determined and optically determined Earth orientation parame- 
al. 2009) is defined by the precise positions of selected compact ters. The problem with these solutions is that since the indirect 
extragalactic radio sources (active galactic nuclei, AGNs) regu- links use stellar sources for the connection between the radio 
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and the optical reference frames, the quality of the link degrades 
with time due to the uncertainties in the measured stellar proper 
motions. 

However, with the sensitive next-generation space astro me- 
rry mission, the European Space Agency's Gaia spacecraft (e.g. 
Mignard 2002; McCaughrean 2012) to be launched in the sec- 
ond half of 2013, a quasi-inertial reference frame can directly be 
established by around 2020, based on measurements of a large 
number of extragalactic sources (probably tens of thousands of 
primary objects) in the optical as well. Gaia will detect a total 
of ~500 000 quasars brighter than the limiting magnitude of 20, 
with a precision similar to ICRF2 (<0.1 milliarcseconds, mas, 
Chariot & Bourda 2012). It will become possible to directly link 
the radio and optical reference frames using a large number of 
common objects for the first time. For consistency, it is impor- 
tant to make the alignment with the highest possible accuracy, 
which requires common objects with excellent optical and ra- 
dio astrometric properties (e.g. Fey et al. 2001). As detailed by 
Bourda et al. (2008, 2010, 201 1), additional observations of new 
radio sources and the construction of the ICRF3 are needed for 
the best possible link, since the potentially most suitable Gaia 
sources for the alignment will not necessarily be the best ICRF2 
sources. The effect of secular aberration drift, i.e. the apparent 
proper motion of extragalactic sources caused by the rotation of 
the solar system around the Galactic centre (Kovalevsky 2003; 
Titov et al. 2011), also has to be taken into account when con- 
structing the next ICRF (Liu et al. 2012). 

The direct connection between the Gaia celestial reference 
frame and the ICRF is essential not only for astrometry, but for 
astrophysical reasons as well. It will become possible to accu- 
rately study the coincidence between the radio and optical emis- 
sion peaks of AGNs on a sub-mas scale (e.g. core shift, Lobanov 
1998; Kovalev et al. 2008), and the observations of the AGNs 
at different wavelengths will only be interpreted correctly if the 
measurements are expressed in a consistent system (e.g. Ivezic 
et al. 2002; Kimball et al. 2010). 

While waiting for the Gaia extragalactic reference frame to 
be constructed, we can perform case studies for the direct ref- 
erence frame link using currently available large optical sky 
surveys. The SDSS (York et al. 2000) is currently the largest 
sky survey available, and although not an astrometric catalogue, 
the faint limiting magnitude (Va;22mag for 95% completion) 
makes it possible to identify the counterparts of many radio- 
loud AGNs that have accurate radio positions. There have been 
several studies cross-referencing the SDSS optical positions di- 
rectly with the radio positions in the ICRF (Frey et al. 2006a; 
Lambert et al. 2006; Souchay et al. 2008; Damljanovic et al. 
2012), or comparing other radio and optical properties based 
on a common optical-radio sample (Ivezic et al. 2002; Kim- 
ball et al. 2010) using e.g. the Very Large Array (VLA) Faint 
Images of the Radio Sky at Twenty-centimeters (FIRST) survey 
catalogue (Becker et al. 1995). Since the astrometric reductions 
of SDSS can be traced back to the Hipparcos catalogue, it can 
also be used to locally connect the radio and optical reference 
frames directly by finding common sources in the SDSS and ra- 
dio catalogues. The ICRF positions of these optical-radio pairs 
are then determined precisely through relative VLBI astrometry 
using nearby ICRF sources as phase-reference calibrators (e.g. 
Frey et al. 2008). In addition, there are other optical catalogues 
that are sensitive enough to make direct connection of quasar po- 
sitions possible between the radio and the optical (e.g. Assafin et 
al. 2005; Souchay et al. 2008), such as the US Naval Observa- 
tory (USNO) CCD Astrograph Catalogue (UCAC3, Zacharias 
et al. 2010), the 2-degree Field Quasar Redshift Survey (2QZ, 



Croom et al. 2004), or other quasar catalogues compiled using 
the above-mentioned radio and optical surveys, e.g. the Large 
Quasar Astrometric Catalogue (LQAC-2, Souchay et al. 2012) 
or the Veron-Cetty & Veron (2010) catalogue. Finally, these 
studies not only provide an independent assessment of the astro- 
metric accuracy of the optical surveys (e.g. Souchay et al. 2008), 
but the comparison of the optical and radio positions could re- 
veal some new information about the physical properties of these 
AGNs. 

In this paper, we present a case study of a direct astrometic 
link between the radio and the optical bands by comparing the 
VLBI positions of AGNs in the ICRF2 with those in the 9 th Data 
Release (DR9) of the SDSS (Ahn et al. 2012). The astromet- 
ric properties of the two catalogues are reviewed and our sample 
selection process is explained in Sect. 2. Section 3 gives the 
statistical characterization of the optical-radio coordinate differ- 
ences. We find a large number of sources that are significantly 
(>3cr) offset between the radio and the optical. These sources 
are presented in Sect. 4. Possible reasons behind the "outliers" 
are discussed in Sect. 5, including astrometric errors and various 
potential astrophysical explanations. Conclusions are drawn in 
Sect. 6. 



2. Selecting common sources in ICRF2 and SDSS 
DR9 

2.1. Astrometric properties of the SDSS Data Releases 

The SDSS DR9 catalogue 1 (Ahn et al. 2012) is the newest re- 
lease of the SDSS-III campaign (Eisenstein et al. 201 1), which is 
an extension of the previous SDSS-I and SDSS-II projects (York 
et al. 2000). It covers ~ 14 500 square degrees of the sky in the 
optical, mainly in the northern galactic hemisphere between right 
ascension 7 h < a <18 h and declination -5°< 6 <+70°. It also cov- 
ers a smaller region in the southern galactic hemisphere between 
21 h < a <4 h and -15°< 6 <+35°. The large sky coverage and the 
faint (V~22 mag) limiting magnitude make it possible to identify 
the counterparts of many radio-loud AGNs that have accurate ra- 
dio positions available in the ICRF2. The astrometric calibration 
of SDSS is described in detail by Pier et al. (2003). The source 
positions derived from the r photometric CCDs (/l r =6165A) are 
calibrated using the 2 nd release of the USNO CCD Astrograph 
Catalogue (UCAC2, Zacharias et al. 2004) and the UCAC rl4 
catalogue (a supplemental set of UCAC at declinations above 
41°). Proper motions are derived from the SDSS+USNO-B cat- 
alogue (Munn et al. 2004). There are ~2-3 magnitudes of over- 
lap between UCAC and unsaturated stars on the r photometric 
CCDs. The UCAC observations are based on the Tycho-2 refer- 
ence stars (H0g et al. 2000), which are directly connected to the 
Hipparcos Reference Frame (HRF). Pier et al. (2003) declares 
that the global astrometric precision of SDSS relative to HRF 
is 45 mas rms per coordinate at r»20 mag and approximately 
100 mas rms at ra;22mag, with additional systematic errors of 
less than 20 mas resulting from the reference catalogues. 

Since DR9 is currently the newest, the most accurate, and 
has the largest coverage among the SDSS data releases, we used 
it for cross-referencing with the VLBI-measured reference po- 
sitions in ICRF2. It also includes the first-year data of the new 
Baryon Oscillation Spectroscopic Survey (BOSS) spectographs 
(Dawson et al. 2013), which focus on obtaining new spectra 
of galaxies in the redshift range 0.15<z<0.8 and quasars with 
2.15<z<3.5, and have a much denser and larger coverage than 

1 http://www.sdss3.org/dr9/ 
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Fig. 1. Graph of the number of counterparts to ICRF2 sources found 
in SDSS DR9 as a function of the search radius. The continuous curve 
represents the total number of hits (AO. whereas the dot-dashed curve 
shows the results corrected with the false match ratio derived from 
Monte Carlo simulation (N—pfN) detailed in Sect. 2.3. The growing 
deviation of the two curves is caused by the increase of false identifica- 
tions at larger search radii. The curves nearly coincide up to 500 mas 
search radius, the value used in this case study, therefore the number of 
false optical-radio identifications is negligible. 



the previous spectroscopic surveys of SDSS-I and SDSS-II de- 
scribed by Abazajian et al. (2009). 



2.2. Astrometry of ICRF2 

The fundamental celestial reference system adopted by the IAU 
has been based on the radio coordinates of AGNs since 1998 
(Feissel & Mignard 1998). It was first realized by ICRF1 (Ma 
et al. 1998) based on the radio positions of 212 extragalactic 
sources distributed over the entire sky. Due to the accumu- 
lating observational data and significant developments and im- 
provements in astrometric VLBI sensitivity and quality, ICRF 
was redefined by the International Earth Rotation and Reference 
Systems Service (IERS) and the International VLBI Service for 
Geodesy and Astrometry (IVS) in 2009. The resulting ICRF2 2 
is currently the realization of the celestial reference system at 
the radio frequencies and is described in detail by Fey et al. 
(2009). It contains the precise positions of 3414 compact radio 
sources, measured with the VLBI technique, with a noise floor 
of only 40 microarcseconds (pas) and an axis stability of 10 pas. 
The coordinate system is maintained using a set of 295 defining 
sources selected on the basis of positional stability and the lack 
of extensive intrinsic source structure. Of the total number of 
sources, 2197 were observed only in the Very Long Baseline Ar- 
ray (VLBA) Calibrator Survey (VCS, Petrov et al. 2008, and ref- 
erences therein), most of them in only one VCS session. These 
sources are located at 5>-30° and their coordinates in most cases 
are as accurate as those of the non- VCS sources, which cover the 
whole sky uniformly. The precision of the source coordinates in 
ICRF2 is better than 1 mas in both right ascension and declina- 
tion, ranging between ~0. 1-0.5 mas depending on the different 
surveys and the number of observations (Fey et al. 2009). 
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Fig. 2. Sky plot (equatorial coordinates in Aitoff projection) of the 
1297 AGNs found in the SDSS DR9 based on a 500-mas radius search 
using the ICRF2 radio source catalogue as reference. This is the sample 
used in our case study of the direct link between the radio and opti- 
cal reference frames. The objects are distributed uniformly in the area 
covered by SDSS DR9. 



2.3. Source selection 

Since the astrometric precision of ICRF2 is more than two or- 
ders of magnitude better than that of the SDSS, we considered 
the radio coordinates as the accurate positions of the AGNs 
when selecting their optical counterparts from the SDSS cata- 
logues. This selection was accomplished by defining a search 
radius around the supposedly "error-free" radio positions and 
identifying the optical sources inside these circles. This selec- 
tion method has been used several times in the past with various 
cut-off radii, between 300-1000 mas, when comparing ICRF- 
SDSS positions (e.g. Frey et al. 2006a,b; Lambert et al. 2006; 
Souchay et al. 2008) and 1'.'5— 3" when using the FIRST ra- 
dio catalogue (e.g. Ivezic et al. 2002). The value chosen for 
the search radius matters because a smaller-than-ideal radius ex- 
cludes sources that would otherwise have optical counterparts. 
A larger-than-ideal radius would contaminate the sample with 
chance identifications, i.e. sources that are in close proximity to 
the AGNs on the celestial sphere, but have no physical relation 
to them. 

As a starting point, the size of the cut-off should be at least 
as big as the ~2-3<x astrometric precision of the less accu- 
rate catalogue, in our case ~200 mas. To determine the ideal 
search radius, we performed a Monte Carlo simulation to cal- 
culate the probability of chance coincidences (i.e. false iden- 
tifications) as a function of the radius used. We constructed 
eight false radio source lists by simply shifting a and 6 for all 
ICRF2 sources by large arbitrary amounts, +(1,2,3,4)°. We 
then tried to find SDSS optical counterparts for these fake "ob- 
jects" (over 27 000 in total) using various search radii ranging 
from 200-3000 mas, with 100 mas increments. The simula- 
tion showed that the probability of chance coincidences (pf) in 
DR9 is -0.01% for 200 mas, remains <0.1% up to 600 mas, 
and reaches ~3% at 3000 mas. The simulations for SDSS DR7 
and DR8 produced similar results and generally agree well with 
Ivezic et al. (2002) and Souchay et al. (2008), who compared 
other catalogues. To find the cut-off radius to use for this case 
study, we also searched for optical counterparts around the real 
ICRF2 positions, and determined the N number of matches as a 
function of the radius. Using the previously obtained pf proba- 
bilities, we subtracted from N the number of false matches, pfN. 
This provided a count corrected for the contaminating false iden- 
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Fig. 3. Normal probability plots of the right ascension (left) and declination (right) optical-radio coordinate differences. The data are normally 
distributed because the graphs approximately coincide with the 45° reference line. There are considerably more objects at both ends of the ordered 
sample data, however, which indicates a heavy tail in the distribution. The curves surrounding the quantile points indicate the 95% confidence 
interval. The quantiles of the sample data are normalized to the same length reference distribution of N(0, 1). 



tifications. Figure 1 shows the number of total hits (AO and the 
corrected hit count (N-pfN) as a function of the search radius. 
This indicates that the ideal cut-off radius is around 500 mas, 
since the sample of the optical-radio matches is the largest here 
without many chance coincidences. 

Among the all-sky set of 3414 ICRF2 sources, optical coun- 
terparts of 1297 (-38%) were found in the SDSS DR9 within 
the search radius of 500 mas (Fig. 2). This is consistent with 
the ~35% sky coverage of the optical catalogue and indicates 
that practically all of the ICRF2 radio AGNs do have an optical 
counterpart with the SDSS limiting magnitude. The identified 
objects are evenly dispersed within the region covered by DR9 
(see Sect. 2.1). In the sample data set, the probability of chance 
coincidences is only ~0.06%, i.e. less than 1 false optical identi- 
fication is expected. This provides us with a clean optical-radio 
AGN sample since these optical counterparts are real identifica- 
tions. Among the 1297 sources, 233 are classified as extended 
(i.e. galaxies) and 1064 as point-like (i.e. quasars) in the SDSS 
DR9. All are primary objects, i.e. their position is from the best 
run in case of multiple observations. The optical coordinates 
are derived from the r photometric CCDs, their average appar- 
ent magnitude is r»;18.9mag. Using this sample of common 
sources in the ICRF2 and SDSS DR9 catalogues, we determine 
and analyse the offsets between the optical and radio positions 
in Sect. 3. 

3. Statistical characterization of optical-radio 
positional differences 

Using the sample data set of the 1297 optical-radio (a, 6) coor- 
dinate pairs, we calculated the angular differences between the 
SDSS DR9 and ICRF2 positions. As mentioned in Sect. 2.2, 
the radio positions can be taken as the reference for the AGNs, 
thus the optical [1] minus radio [2] coordinate differences can be 
calculated using spherical trigonometry formulae as 

Aa = (a i - a%) cos 62 , (1) 



AS = Si - 6 2 , (2) 

cos A = cos(ai - (X2) cos 61 cos 62 + sin 6\ sin 62 , (3) 

where Aa and Ad are the differences in right ascension and decli- 
nation, respectively, and A is the total angular difference between 
the optical and radio positions. 

3. 1 . Testing of normality 

To characterize statistically the optical-radio positional dif- 
ferences, we tested the normality of the calculated random 
variables. This was accomplished by constructing a normal 
quantile-quantile (Q-Q) plot for both Aa and A6 (Fig. 3). Nor- 
mal Q-Q plots are probability plots that compare a sample data 
of unknown distribution with the standard normal distribution 
N(Q, 1). A Q-Q plot is commonly used for comparing a data set 
to a theoretical model and is a robust graphical method of distri- 
bution analysis (e.g. Fisher 1983; Evans et al. 2000; Rosenkrantz 
2000; Das & Resnick 2008). It is also used in astronomical data 
interpretation (e.g. Konig & Timmer 1997; Pestana & Cabrera 
2004; Huff & Stahler 2008). In our study, the general trend of the 
points is somewhat flatter in both cases than the 45° line, which 
indicates that the middle part of the data is slightly less dispersed 
than an N{0, 1) distribution. Moreover, as seen from the arcs at 
both ends on the probability plots (Fig. 3), the distributions of 
our data have heavier tails than a normal distribution. However, 
since the quantile pairs on both normal Q-Q plots approximately 
lie along the line y — x, both the right ascension and declination 
differences can be considered to have a normal distribution with 
a fat tail. The correlation coefficient between the two quantities 
is very low, 0.017. This means that Aa and A6 can, for all things 
considered, be treated as independent normal random variables. 
Because the involved angular distances are very small, the total 
positional offset can be approximated with planar trigonometry 
as A« VAa 2 + Ad 2 , so we can assume a Rayleigh distribution 
for its statistical characterization. 
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Fig. 4. Histogram of optical- VLBI total positional differences for 
1297 ICRF2 sources identified in the SDSS DR9 within a search radius 
of 500 mas. The distribution is consistent with ~57 mas (lcr), assuming 
Rayleigh distribution. The vertical dashed line indicates the 3<r value 
of 170 mas. A total of 51 sources have optical-radio offsets larger than 
3cr. 



3.2. Data weighting and statistical parameters 

When determining the Gaussian parameters of Aa and AS, it is 
possible to refine the calculations by introducing weights charac- 
terizing the reliability of the optical positions. SDSS uses algo- 
rithms (described in Petrosian 1976; Blanton et al. 2001; Yasuda 
et al. 2001) to automatically distinguish extended sources (i.e. 
galaxies) from point-like objects (i.e. stars or quasars). However, 
as detailed in Ahn et al. (2012), some inconsistencies in classifi- 
cation have been found in earlier SDSS Data Releases. To avoid 
any possible problem, we did not rely on the galaxy/quasar clas- 
sification of SDSS when calculating the statistical characteristics 
of the offsets between the radio and optical positions. Instead, 
we took the Petrosian radius (Blanton et al. 2001; Yasuda et al. 
2001), a measure of the optical size of a given source, and used 
it as a reciprocal weight for the individual source coordinates. 
This way we could take into account that point-like sources have 
more reliable position measurements. At the same time, we did 
not exclude extended optical sources as in an earlier study with 
ICRF2 and SDSS DR7 (Orosz & Frey 2012). 

The distribution of the A positional differences of SDSS- 
ICRF2 AGNs is shown in the histogram of Fig. 4. It is gener- 
ally consistent with a 57-mas positional uncertainty, very close 
to the declared global SDSS astrometric precision (Pier et al. 
2003). The weighted optical minus radio right ascension differ- 
ences (Aa) have <x Aq .=55 mas standard deviation and a negligible 
/UAa= _ 3 mas mean value. The weighted declination differences 
(AS) have similar standard deviation of <x Aa =54 mas and also a 
negligible fi^ = 4 mas mean value. As found in the normality 
test, the distribution has a fat tail, which means that there are sig- 
nificantly more sources with >3cr positional offset than expected 
statistically. The theoretical distribution would only give ~four 
sources (0.27%) above the 3cr level, as opposed to the detected 
51 outliers, i.e. AGNs with optical-radio separations exceeding 
170 mas. They represent ~4% of the sample. Furthermore, ~1% 
of the sources are even beyond the 5cr mark, which should be 
practically impossible in a population with standard normal dis- 
tribution. Earlier studies using SDSS DR4 data with 524 match- 
ing quasars and galaxies (Frey et al. 2006a), SDSS DR5 data 
with 735 matching objects (Frey et al. 2006b), and SDSS DR7 



data with 806 matching quasars (Orosz & Frey 2012) led to sim- 
ilar results, with outlier ratios also around 4%. 



4. Sample of positional outlier objects 

Comparing the radio positions of ICRF2 sources to their opti- 
cal counterparts in SDSS DR9, we found a total of 51 AGNs 
where the optical-radio offset exceeds the 3<x level determined 
in Sect. 3. These sources apparently show an even distribu- 
tion in the sky within the coverage of SDSS DR9. Except for 
J1506+3730, none of them is a defining object of ICRF2. The 
complete list of positional outliers is presented in Table 1 . For 
each AGN, Col. 1 gives a short IAU designation constructed 
from the J2000.0 equatorial coordinates as JHHMM+DDMM. 
Note that the format of the complete ICRF2 designations that 
can be derived from the right ascension and declination is ICRF 
JHHMMSS.s+DDMMSS (Fey et al. 2009). Based on informa- 
tion from the literature, in Sect. 5.5 we comment on the prop- 
erties of seven objects marked here with asterisks. The J2000.0 
equatorial coordinates are taken from the ICRF2 catalogue (Fey 
et al. 2009) and listed in Cols. 2 and 3. The total optical-radio 
angular offsets (A), and the offsets broken down to coordinate 
components (Aa and AS) are given in Cols. 4, 5, and 6, respec- 
tively. The apparent r magnitudes (m T ) taken from SDSS DR9 
can be found in Col. 7. Redshifts (where available) and a simple 
optical classification from SDSS DR9 (Q: quasar, G: galaxy) are 
given in Cols. 8 and 9. Column 10 provides the references for 
redshifts. 

We found spectroscopic redshift measurements in the liter- 
ature for 27 out of 51 sources. These redshifts range between 
0.04^z^3; half of them are below 0.6, and about a quarter of the 
sources have z^2. We examined whether there is a connection 
between the angular separation of the optical-radio positions and 
the redshifts, and found no evidence for it. About one third of 
the outliers are classified as quasars in SDSS DR9, with the rest 
being galaxies, which generally have lower redshifts in the sam- 
ple. Remarkably, there are three known cases of gravitational 
lenses in our list (J0134-0931, J0414+0534, and J1601+4316). 
For these sources, the redshifts of both the lensing objects and 
the lensed background sources are given in Table 1 . Two other 
AGNs (J 1006+3454 and J 1301 +4634) may show double-peaked 
narrow [O III] emission lines among the total of ten sources that 
have optical spectra available in the SDSS DR9. There are also 
a couple of objects with large extended radio structures or other 
interesting properties, which may well be related to the signif- 
icant optical-radio positional differences. In what follows, we 
briefly review possible explanations of the large positional off- 
sets, and discuss the properties of the most interesting individual 
sources in our sample. 



5. Discussion 

When comparing directly the radio and optical positions of 
AGNs, we naturally assume that the optical and radio emission 
peaks physically and spatially coincide. Because the activity of 
these distant extragalactic sources is driven by matter accretion 
onto their central supermassive black holes and is confined to 
their close vicinity, this seems a plausible first approximation. 
Theoretically, the apparent origin of the inner radio jet, usually 
called the VLBI core, depends on the observing frequency due to 
opacity effects (see e.g. Lobanov 1998). Actual VLBI measure- 
ments of this core shift indicate that it can be as high as ~ 1 mas 
for certain sources (e.g. Kovalev et al. 2008; Sokolovsky et al. 
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2011) . Kovalev et al. (2008) estimate the average shift between 
the cm-wavelength radio core and the optical core as ~0.1 mas. 
This is more than two orders of magnitude lower than the SDSS 
positional accuracy. The core shift is therefore negligible and 
cannot be the cause of large optical-radio outliers in our sample. 
However, this effect should be considered for the alignment of 
the future Gaia optical reference frame with the ICRF. Similar 
arguments are valid for possible optical photocentric variability 
of quasars (Popovic et al. 2012): while the effect can be as large 
as several mas for low-redshift AGNs, it is certainly negligible 
in our case. In the next subsections, we look into some causes 
that could possibly be behind the positional outliers we found. 

5.1. Errors in positions 

Before discussing possible astrophysical explanations, we inves- 
tigate whether systematic astrometric calibration errors in the 
SDSS database can cause the occasionally large optical-radio 
positional offsets. In the course of our work, we performed sim- 
ilar analyses to find optical counterparts to ICRF2 radio sources 
by cross-referencing their positions with the SDSS DR7 (Abaza- 
jian et al. 2009) and DR8 (Aihara et al. 2011a), and calculated 
the optical-radio coordinate differences. The SDSS DR7 cata- 
logue covers ~ 12 000 square degrees in the sky, mainly in the 
northern Galactic hemisphere and around the equator. Its astro- 
metric calibration is similar and the accuracy is basically iden- 
tical with DR9 (Pier et al. 2003), apart from some proper mo- 
tion errors for stars at low Galactic latitudes later corrected in 
DR9. Among the 51 positional outlier sources we found in DR9 
(Table 1), 37 are located and identified as photometric objects 
in DR7 as well, reflecting its smaller sky coverage. We used 
this overlapping subsample to examine the consistency of out- 
liers between the two databases, and found that 27 of the 37 
sources are outliers in both DR7 and DR9. None of them is 
associated with the reportedly miscalibrated runs in DR7 (Ahn 
et al. 2012). This 73% overlap between the two lists suggests 
that although the majority of our outliers listed in Table 1 are 
robust detections, there might be cases where small (local) as- 
trometric calibration issues in SDSS contribute to the positional 
offsets. Worth noting is that the optical counterparts of ICRF2 
radio sources are typically faint (r=18.88±1.83 mag), close to 
the SDSS limiting magnitude, therefore their optical positions 
are less precisely determined. These problems will be alleviated 
when the more accurate Gaia catalogue becomes available for 
such an optical-radio study. 

We repeated our analysis with DR8, which has the same sky 
coverage as DR9. However, when calculating the offsets be- 
tween the optical and radio coordinates, we found an anomaly 
in the declination differences. Almost all sources above i5~40° 
consistently showed an offset of ~260 mas relative to the ICRF2 
coordinates. This indicated a systematic error in the astromet- 
ric calibration of SDSS DR8. Our independent result is consis- 
tent with what is reported and detailed by Aihara et al. (201 lb). 
This astrometric calibration problem has been fixed in the SDSS, 
and also prompted more rigorous astrometric quality-assurance 
measures in DR9 using a set of reference catalogues (Ahn et al. 

2012) . This justifies the use of DR9 for our study. According 
to our normal probability plots in Fig. 3, there are no dominant 
systematic components in the optical-radio offsets using DR9. 
(The same is valid for DR7.) This means that there are no global 
errors in the astrometric reductions of these data releases . The 
even distribution of the outlier AGNs in the DR9 sky coverage 
(Fig. 2) also supports this notion. 



5.2. Errors in identifications 

As we have shown in Sect. 2.3, our source selection method with 
a 500 mas search radius around the positions of ICRF2 objects 
guarantees that the matched optical-radio AGN sample is prac- 
tically free from false identifications, i.e. physically unrelated 
radio and optical sources. This is because the average distance 
between quasars in the SDSS is much larger, in the order of ar- 
cminutes (cf. Palanque-Delabrouille et al. 2013). However, it is 
in principle not excluded that the radio position of a particular 
object does not refer to the AGN core (which is in fact the base 
of the inner jet, close to the central engine) but to a brighter and 
compact component in one of the outward moving relativistic 
radio jets. This is certainly not the case for the ICRF2 sources 
in general, since these are among the most prominent compact 
radio AGNs, often with sensitive multi-frequency radio imag- 
ing observations. But for at least one of our outlier sources, a 
peculiar quasar with a complex extended two-sided jet structure 
(J1526-1351), the ICRF2 catalogue indeed contains the position 
of a jet component is actually brighter than the quasar core (see 
Sect. 5.5 for the details). 



5.3. Gravitational lensing 

There could be astrophysical explanations for optical-radio off- 
sets of ~ 100 mas as well. Strong lensing caused by the interven- 
ing gravitational potential of a foreground galaxy can produce 
multiple and distorted images of a background object. The sepa- 
rations of gravitationally lensed AGN images are typically in the 
order of 1" or less. Thus the optical images are blended, possi- 
bly together with the lensing galaxy, and this complex structure 
remains unresolved in the SDSS. The optical position refers to 
the photocentre. On the other hand, the radio coordinates de- 
rived from VLBI data with much higher angular resolution usu- 
ally refer to a particular (the brightest and most compact) image 
of the gravitationally lensed source. A special scenario where 
the lensed quasar is radio-loud but optically faint while the lens- 
ing galaxy is optically bright cannot be excluded. This could 
also lead to sub-arcsec apparent separation between the optical 
and radio positions. 

Successful systematic searches for gravitational lenses in the 
SDSS database were performed using spectroscopic data (e.g. 
Bolton et al. 2008) or morphological and colour-selection crite- 
ria (e.g. Inada et al. 2012). The candidates were then followed- 
up by higher-resolution imaging and spectroscopic observations 
for verification. Recent simulations show that Gaia will detect 
~0.6% of the quasars that are multiply imaged due to gravi- 
tational lensing (Finet et al. 2012). Here we deal with radio- 
selected AGNs taken from the ICRF2 catalogue. To estimate 
the probability of finding lensed objects in our case, we there- 
fore applied the result from the largest sample of strongly lensed 
flat-spectrum radio sources in the Cosmic Lens All-Sky Survey 
(CLASS, Browne et al. 2003), where -0.14% was found for the 
point-source lensing rate. In sharp contrast with the statistical 
expectation, 3 of our 51 outlier sources (~6%) are known grav- 
itational lens systems (see Sect. 5.5 for details). This indicates 
that the method of finding AGNs with significant optical-radio 
positional offsets may be an efficient tool for identifying gravi- 
tational lens candidates. Note that for many of the objects in Ta- 
ble 1, optical spectroscopic data are unavailable at present, and 
sub-arcsecond separation gravitational lensing cannot be con- 
vincingly excluded. 
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5.4. Interacting active galactic nuclei 

Another possible astrophysical cause for the measured optical- 
radio positional offsets could be related to dual AGNs. The pres- 
ence of dual accreting supermassive black hole systems is a nat- 
ural consequence of hierarchical structure formation in the Uni- 
verse through mergers of galaxies (e.g. Begelman et al. 1980). 
Hydrodynamical simulations (Van Wassenhove et al. 2012) sug- 
gest that simultaneous AGN activity is mostly expected at the 
late stages of mergers, at or below ~1-10 kpc separations. 
Assuming a cosmological model with H^-IQ km s _1 Mpc -1 , 
Q m =0.3, and Qa=0.7, 1 kpc linear size corresponds to SS300 mas 
angular size at any plausible redshift above z=0.2, or smaller if 
the dual system is inclined to the plane of the sky. Thus the 
expected angular separations are typically below the resolution 
limit of SDSS, but comparable to the optical-radio positional 
offsets found in this paper. For example, a dual AGN system of 
which one component is presently in its short-lived radio-loud 
activity phase but the other one is not, would easily result in a 
detectable positional offset. 

The SDSS spectroscopic database was used to find candidate 
dual AGNs via searching for double-peaked narrow [O III] or 
other emission line profiles (e.g. Smith et al. 2010). These may 
indicate gravitationally bound dual AGNs with distinct narrow- 
line regions (NRL). In this model, the two components have dif- 
ferent radial velocities due to their orbital motion, resulting in 
the doubling of the line profiles. However, peculiar gas kine- 
matics and jet-cloud interaction in a single NLR can also lead 
to similar spectral signatures. Indeed, confirmed kpc-scale dual 
AGNs seem to add up only a few percent of the candidates with 
double-peaked [O III] emission lines (Shen et al. 201 1). Accord- 
ing to the estimates of Rosario et al. (201 1) and Fu et al. (201 1), 
only ~0.3% or less of the low-redshift (z^0.6) SDSS quasars 
host dual accreting black holes separated on kpc scales. This 
ratio decreases with increasing redshifts (Yu et al. 201 1). 

A sample of objects with large optical-radio positional off- 
sets and double-peaked [O III] optical emission lines may pro- 
vide better candidates for actual dual AGNs. In our case, only 
10 out of 51 sources have optical spectra available in the SDSS 
DR9. By visual inspection, two of them show indication of dou- 
ble peaks in their [O III] emission lines (see Sect. 5.5), albeit 
with a low signal-to-noise ratio. With more complete spectral 
data, more precise Gaia astrometry in the future, and additional 
observational verifications (e.g. with adaptive optics imaging or 
optical slit spectroscopy, Shen et al. 2011; Rosario et al. 2011), 
this possibility is worth investigating in more detail. 

After final coalescence, a supermassive black hole may be 
kicked out of the centre of its host galaxy. Numerical simula- 
tions of recoiling black holes predict up to kpc-scale separations 
(for a recent review, see Komossa 2012). Observationally, these 
objects could also appear as radio AGNs offset from their host 
galaxies. 

5.5. Comments on individual sources 

Properties that may be relevant for this study are collected from 
the literature for the following individual objects in Table 1. 

J0134-0931. This object is the quintuple quasar, a gravita- 
tionally lensed compact radio source at z s =2.216 consisting of 
five components with a maximum separation of ~0'/7 (Winn et 
al. 2002; Gregg et al. 2002; Hall et al. 2002; Winn et al. 2003; 
Keeton & Winn 2003). The VLBI coordinates refer to one of the 
lensed images, component A. On the other hand, all five images 
of the background quasar, as well as the lensing pair of galaxies 



(zi =0.765), are unresolved in SDSS. It is therefore not surprising 
that we found an offset of 175 mas between the radio position of 
component A and the SDSS optical position, which should refer 
to the peak of the blended emission of the lensing galaxy pair 
and the gravitationally lensed components. 

J0414+0534. This radio source at z s =2.639 (Lawrence et 
al. 1995), also known as 4C +05.19, is gravitationally lensed 
(Hewitt et al. 1992) by an elliptical galaxy (Schechter & Moore 
1993) at zi=0.958 (Tonry & Kochanek 1999). The system of four 
strong lensed radio components separated by nearly 2" has been 
the subject of several VLBI studies (e.g. Patnaik & Porcas 1996; 
Trotter et al. 2000; Ros et al. 2000; Volino et al. 2010). A par- 
tial Einstein ring connecting the three brightest images is seen 
in a deep high-resolution Hubble Space Telescope (HST) image 
(Falco et al. 1997). The ICRF2 coordinates refer to component 
Al, while the optical position is affected by the complex struc- 
ture of the different quasar images and the lens, only partially 
resolved in SDSS. 

J1006+3454. Also known as 3C 236, this prominent radio 
source at z=0.1 is among the largest radio galaxies observed (e.g. 
Willis et al. 1974; Barthel et al. 1985; Schilizzi et al. 2001; Labi- 
ano et al. 2013, and references therein). The overall radio struc- 
ture of this extensively studied object has a narrow morphology 
of ~40' angular size, extending in about SW-NE position an- 
gle. Radio interferometric observations reveal a double-double 
morphology, with a 2-kpc scale compact steep-spectrum (CSS) 
source inside. The accurate VLBI position refers to the compo- 
nent B2 identified with the radio core by Schilizzi et al. (2001). 
The SDSS optical position is offset by 217 mas, roughly along 
the position angle of the radio structure. The optical-radio offset 
is consistent with Fig. 8 of O'Dea et al. (2001) that displays an 
overlay of VLBI and HST images. The apparently renewed radio 
activity in 3C 236, the presence of young star-forming regions in 
the dust lane (O'Dea et al. 2001), and the disturbed kpc-scale op- 
tical morphology suggest a galaxy merger that might have led to 
the reignition of the radio AGN (see e.g. Labiano et al. 2013, 
and references therein). Notably, the SDSS DR9 spectrum of 
J1006+3454 (3C 236) hints at a double-peaked narrow 5007 A 
[O III] emission line. 

J1301+4634. A low-redshift (z=0.2) galaxy, which is the 
central, brightest member of a small cluster selected from SDSS 
photometric data (Koester et al. 2007). Visual inspection of the 
SDSS DR9 spectrum suggests broadened or double-peaked nar- 
row [OIII] emission lines. 

J1506+3730. Although the radio-loud AGNs are typically 
found in ellipticals, the host of J1506+3730 is an inclined disk 
galaxy, with an optically obscured nucleus. Significant neutral 
hydrogen and molecular absorption is detected towards the ra- 
dio AGN, arising from a fast gas outflow, a possible result of 
jet-cloud interaction (e.g. Carilli et al. 1997; Kanekar & Chen- 
galur2008, and references therein). The optical-radio positional 
difference we found can be reconciled with the fact that the nu- 
cleus of this red quasar is heavily obscured in the optical. The 
SDSS DR9 optical position is offset from the VLBI position by 
316 mas in the ~SW direction, broadly coinciding with the po- 
sition angle of the radio jet (e.g. Polatidis et al. 1995). 

J1526-1351. This is a high-luminosity CSS quasar, with a 
complex radio structure spanning ~1". Based on polarization- 
sensitive and dual-frequency VLBI imaging observations, Man- 
tovani et al. (2002a,b) found a peculiar non-collinear two-sided 
jet in this source. The flat-spectrum core is ~100 mas south of 
the radio peak, whose position is listed in the ICRF2 catalogue. 
The peak in fact coincides with a bright and compact component 
at the end of the northern jet. It is therefore reasonable to assume 
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that the optical AGN position is closer to the radio core, roughly 
explaining the offset we found. 

J1601+4316. This source is the third case of known strong 
gravitational lenses in our sample, discovered by Jackson et al. 
(1995). The lens is a spiral galaxy at zi=0.414, the background 
source is a doubly imaged quasar at z s =1.589 (Fassnacht & Co- 
hen 1998). The ICRF2 radio catalogue lists the coordinates of 
both the brighter component A and the weaker component B, 
separated by l'/4 (e.g. Koopmans et al. 2000; Patnaik & Kem- 
ball 2001). Our selection algorithm found an optical counterpart 
to component B in SDSS DR9, the one seen very close to the 
lens (Jaunsen & Hjorth 1997). The position angle of our optical- 
radio positional offset (322 mas towards ~NW) is consistent with 
the direction of the lensing galaxy and component A. 

6. Conclusions and outlook 

We performed a case study for directly linking the radio and opti- 
cal reference frames using common objects in the ICRF2 (Fey et 
al. 2009) and SDSS DR9 (Ahn et al. 2012) catalogues. We found 
optical counterparts for 1297 radio-loud AGNs, practically all 
that are located in the SDSS DR9 footprint in the sky, within a 
search radius of 0'.'5 around the accurate VLBI-determined po- 
sitions. Our overlapping ICRF2-SDSS DR9 sample is free from 
false identifications. The optical-radio coordinate differences in 
right ascension and declination follow normal distributions with 
fat tails. This allows us to characterize the overall astrometric ac- 
curacy of SDSS DR9: both equatorial coordinates of the match- 
ing extragalactic objects are in general consistent with no offset 
with respect to the ICRF2, with ~55 mas standard deviation. 

We identified 51 AGNs (~4% of the sample) for which the 
optical-radio positional offset is significant, exceeding 170 mas 
(~3<x). We argued that there is an astrophysical cause behind the 
majority of these outliers. The presence of significant positional 
outliers underlines the importance of a careful reference source 
selection for the precise alignment of the most accurate radio 
reference frame with the future optical frame to be constructed 
from the measurements of the Gaia astrometric space telescope. 

The question of optical-radio positional differences was also 
investigated by Camargo et al. (201 1) using a different approach, 
with a smaller sample. They performed accurate relative op- 
tical astrometric measurements in the fields around 22 ICRF2 
sources. Four of their sources showed offsets relative to their 
ICRF2 positions larger than their 3cr confidence level, ~80 mas. 
Camargo et al. (201 1) concluded that these separations cannot be 
explained merely by statistical fluctuations or systematic errors 
in the optical reference frame, and they might be related to the 
relatively more complex VLBI structure of the given quasars. 

We offered some viable explanations for the positional out- 
liers listed in Table 1. As shown for the peculiar quasar 
J1526-1351 as an example, the ICRF2 position may refer to 
a bright component farther along the radio jet, and not to the 
true VLBI core, which is supposed to be close to the location of 
the central supermassive black hole. Although this situation is 
believed to be rare if not unique, firm proof should come from 
sensitive, high-resolution multi-frequency VLBI observations of 
individual radio reference frame objects, not just typical snap- 
shot imaging. 

There are three known cases of strong gravitational lensing 
in our sample, nearly 50 times more than expected on statistical 
grounds. For these objects, the ~100-mas scale offsets between 
the radio position (which refers to one of the lensed images) and 
the optical photocentre is not surprising. It is possible that there 
are more gravitationally lensed objects in our sample, waiting 



for identification by follow-up spectrosopic observations or sen- 
sitive high-resolution imaging. 

Positional non-coincidences between the optical and the ra- 
dio can also be caused by interacting AGNs containig dual 
AGNs. Observationally, known dual AGNs with ~l-10-kpc 
scale separations are rare, partly because the period of their si- 
multaneous activity within the lifetime of the AGNs is probably 
short. Moreover, it is difficult to confirm their existence due to 
the high-resolution and sensitive imaging and spectroscopic ob- 
servations required. At present, there is no efficient selection 
method to apply for finding firm dual AGN candidates. It re- 
mains to be seen if the significant optical-radio positional differ- 
ence in the case of some radio-loud AGNs found in this paper 
is a good indication of AGN duality. A much larger sample will 
be offered for such a study by the Gaia astrometric catalogue of 
extragalactic sources. Once this becomes available, as Browne 
(2012) suggests, a future extensive radio interfero metric survey 
conducted with the e-MERLIN array in the United Kingdom, 
targeting elliptical galaxies that host radio-loud AGNs, could re- 
veal displacements. These may arise either from dual AGNs or 
recoiling supermassive black holes that are expelled from the 
galaxy after the coalescence of a binary system (Komossa 2012). 
On a short term, sensitive e-MERLIN imaging of our sample in 
Table 1, probing the ~100-mas scale radio stucture, should be 
able to provide clues for tracking down the causes of significant 
positional offsets between SDSS DR9 and ICRF2. Looking for 
signatures of gravitational lensing or double-peaked narrow opti- 
cal emission lines would also require high-quality spectroscopic 
observations of these optically faint sources. 
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Table 1. List of the 51 optical-radio positional outlier objects. 



Short designation 


Position 


in radio 


Offset in optical 


m r 


z 


Type 


Ref. 




RA (J2000.0) 


Dec (J2000.0) 


Total 


RA 


Dec 




(lens, source) 




z 




h m s 


Off/ 


ffltlS 


ffltlS 


ffltlS 


mag 








J0038-2120 


00 38 29.9547 


-21 20 04.023 


225 


117 


192 


18.59 


0.338 


G 


1 


J0041+1339 


aa /i 1 1 t 1 1 1 a 

00 41 17.21 10 


13 39 27.527 


348 


179 


-298 


22.53 




G 




J0106+3402 


01 06 00.2934 


34 02 02.988 


431 


428 


-52 


24.05 


0.579 


G 


2 


|A1 MA_i_9^7Q 

juiuo+zj jy 


M1 MA 1 M QAQM 


9^ 7M /1QA 


353 


-119 


-333 


1 / .41 


M 1 QQ 


/ ■ 

VJ 


9 

z 


10106-031 5 


01 06 43 2287 


-03 1 5 36 296 


207 


44 


-202 


17.65 




n 

V 




j/\ i i yf aai 1 

J0134-0931* 


ai i /i ic ti ti ti ti 

01 34 35.6666 


c\c\ 1 1 r\o oiA 

-09 31 02.879 


175 


36 


-172 


21.35 


(0.765, 2.216) 


G,Q 


3, 4 


J0146+2110 


01 46 58.7839 


r\ 1 1 n /-\ J 101 

21 10 24.381 


277 


122 


-248 


23.03 




Q 




J0216-0118 


02 16 05.6638 


-01 18 03.397 


193 


-162 


-106 


18.73 




Q 




imi/; mim^ 

JUZ10— U1UJ 


AO 1 A 19 9 1 1 Q 

WZ 10 1Z.Z1 iy 


M1 M^ 1 8 89A 
— Ul VJ lo.oZO 


179 


-141 


-111 


1 7 AR 
1 / .Oo 


1 Z1Q9 
1 .+yZ 


ri 

V I 


c 

J 


10334+0800 


03 34 53 3167 


08 00 14 419 

vy vJ vy vy i i . i i ^ 


272 


-179 


-205 


22.49 


1.982 


G 


2 


J0335-0709 


r\ /-> ic C l ~l A C C i"\ 

03 35 57.0552 


r\^i r\r\ c c a e a 

-07 09 55.854 


347 


329 


109 


21.86 




G 




■ /\ At A . a C 1 -1 

J0414+0534* 


04 14 37.7678 


AC 1/1 A r\ 11C 

05 34 42.335 


189 


-189 


-1 


22.72 


/A A C O 1 /^1A\ 

(0.958, 2.639) 


G,Q 


6, 7 


J043 1+2037 


i"v /i ti ai a/" i -i 

04 31 03.7614 


20 37 34.265 


203 


168 


114 


18.04 


0.219 


G 


8 


IMA^j-9^9 
JU+JJ+Zj JZ 


Mzt 3^ 1A ^89Q 

U+ JJ J4.JOZ7 


9« 09 «o AQ7 
ZJ jL jy.ov / 


249 


-238 


75 


9M 9^ 
ZW.ZJ 




V 




10523+6007 


05 23 1 1 0082 


60 07 45 720 

vy V./ W I i^J . 1 i— Vy 


194 


-43 


189 


23.02 




n 

V 




1 A r r- 1 , AH 1 

J0552+0313 


r\c c\ c r\ 1 ai f 

05 52 50.1015 


r\i 1 i n i/ii 

03 13 27.243 


281 


-220 


175 


23.14 




G 




T/\ /I /I . 1 1 

J0644+2911 


/\ f a A A A 01CO 

06 44 44.8158 


i r\ 11 t~\ a f\ 1 o 

29 ll 04.018 


264 


-176 


197 


21.84 




G 




J0729-1320 


AA 1 i—i O 1 

07 29 17.8177 


1 i i r\ r\i ni 

-13 20 02.272 


448 


360 


-266 


16.68 




G 




TM77A_i_9Q^/1 

ju / jO+zyj+ 


M9 7 A 1 7 AA1 1 
U/ JU 1 J. 001 1 


1Q 10 1 
Zy J+ ZZ. loO 


192 


168 


94 


99 AQ 
ZZ.+V 




C\ 

V 




T08 17+3227 

JUU J. / I ~J z— / 


08 1 7 28 5423 

\J\j J. / (J • »y i ~j 


32 27 02 926 


226 


-85 


-210 


21.08 




Q 




IAfl (1 ,1 C 1 

J0843+4537 


08 43 07.0942 


a c n a i Cir\'~i 

45 37 42.897 


177 


163 


-71 


17.57 


0.192 


G 


5 


J0854+6218 


08 54 50.5763 


/'i 1 o c r\ 1 r\ 1 

62 18 50.191 


176 


-115 


-134 


18.26 


0.267 


G 


5 


J0902+4310 


09 02 30.9200 


/II 1A 1/1 1 

43 10 14.166 


228 


-201 


108 


20.57 


2.41 


Q 


9 




1 M MA M1 H^CYX 
1U UO Ul . / JUJ 


7/1 ^/l 1 M /1M1 
J)+ J+ 1U.+U1 


217 


-180 


120 


1 ^ MQ 

i j.uy 


M MQQ 




c 

J 


11022+4239 


10 22 13 1323 


42 39 25 612 

Z— ^ / .y a-j^j • vy i z_ 


260 


259 


-23 


22.59 


0.991 


Q 


2 


ii AOO i OfiTC 

J1033+3935 


10 33 22.0610 


an 1C C1 AOO 

39 35 51.083 


259 


-259 


14 


21.45 


1.095 


G 


10 


Jl 150+4332 


11 50 16.6027 


/i i ii r\c /■ 

43 32 05.906 


202 


183 


-86 


20.94 


3.037 


Q 


1 1 


J1254+0859 


12 54 58.9577 


f\Ci C l\ A i~l C A r\ 

08 59 47.549 


239 


-230 


-65 


23.00 




G 




J 1 Jul ++0 J+* 


11 A] ^9 AMA^ 
LJ Ul JZ.OUOJ) 


ztA 'XA M9 QAM 
+0 J+ WZ.y+W 


202 


130 


-154 


1 A II 
10. J J 


M 9MA 


V i 


c 

J 


T1312+2531 


13 12 14 2889 


25 31 13 175 

A-c^j j i J. ^y . J. / -y 


336 


180 


284 


23.32 




G 




J1312+4828 


i o 1 o icon 

13 12 43.3537 


/lO lO TA Ayl 1 

48 28 30.941 


411 


395 


-112 


20.95 


0.501 


G 


2 


J1313+6735 


13 13 27.9863 


f —1 r~t C C f\ lOI 

67 35 50.382 


439 


-1 


-439 


21.92 




G 




11 ,11 „ /I C C /I 

J1414+4554 


1/1 1/1 1/1 O C l"> 

14 14 14.8526 


AC C A A Ci 

45 54 48.720 


171 


-164 


47 


20.16 


0.186 


G 


11 


11 ztAMo-M1 97 
J 1++U+U1Z / 


1 A ztM Azt7M 
1+ +W JJ).0+/W 


M1 97 M<; 9 1 M 


262 


262 


1 


91 MA 
Z 1 .WO 




n 

V I 




11451+1343 


14 51 31 4910 

J. T ^y _L J 1 . T y J. vy 


13 43 24 001 


420 


234 


349 


21.86 




G 




T1 AC/" . c f\ a a 

J1456+5048 


1 /I tT* O 1 1 f\^1 

14 56 08.1197 


CA ( O I/' 1AA 

50 48 36.300 


228 


212 


-85 


22.87 


0.480 


Q 


12 


11 r oo . aai n 

J1503+0917 


15 03 00.8995 


r\r\ 1 <—i c a r\ a i 

09 17 58.983 


291 


-290 


24 


21.98 




G 




J1506+3730* 


1 c /"\ /i r\r\ ctaa 

15 06 09.5300 


n ta c 1 111 

37 30 51.133 


316 


-237 


-209 


21.23 


0.672 


G 


13 


11 ^9A 1 7^1 .j- 
J 1 jZO— 1 JJ 1 * 


1 ^ 'nQ a Ann 


17^1 MM 1 A/1 
Ij Jl UU. 10+ 


192 


-78 


-176 


1 Q AA 
ly. ++ 


1 AC7 
1 .08 / 


V 


1 A 
1 + 


11 543+0452 

J i ^y r ~J i v/ T^y z^ 


IS 43 33 Q?58 

J. -y i ^> ^> • J a-, \j 


04 52 19 320 

Vy ^7 z— 1 y . ^ ' z— V./ 


413 


409 


-58 


13.60 


0.040 


G 


5 


T 1 / /\ 1 1 T 1 i" 

J1601+4316* 


16 01 40.5154 


/II 1 /< /I r-7 ^7 

43 16 46.477 


322 


-139 


291 


20.84 


/A AAA 1 COA\ 

(0.414, 1.589) 


G,Q 


15 


11 /"A1 1 C C ,1 

J1603+1554 


16 03 38.0619 


ic c ,1 r\i ice 

15 54 02.355 


205 


87 


185 


15.22 


O.llO 


G 


5 


11 /"A < 1 f 

J 1604+ 1926 


16 04 49.9938 


t r\ i /~\/ii 

19 26 20.942 


194 


109 


-160 


22.36 




G 




T1 A0^j./I 1 7/1 
J lOZj+41 J+ 


1 A 9^ ^9 AAQ9 

to zj j / .ooy / 


/l 1 7/1 /1M AOQ 
+ 1 J+ +U.0Z!7 


259 


109 


235 


99 7^ 
ZZ. JJ 


9 ^ 
Z.JJ 


/ ■ 

VJ 


1 A 
14 




1 6 48 01 5356 


22 24 33 148 


179 


-23 


178 


21.70 


823 


r; 


■J 


J2052+1619 


20 52 43.6199 


16 19 48.828 


231 


229 


31 


21.82 




Q 




J2150+1449 


21 50 23.6071 


14 49 47.895 


189 


-38 


185 


21.84 




Q 




J2210+0857 


22 10 06.0503 


08 57 29.564 


193 


-86 


173 


18.84 




G 




J2259-0811 


22 59 00.6888 


-08 11 03.043 


235 


-231 


-41 


20.32 


1.380 


Q 


5 


J2346+3011 


23 46 46.2508 


30 11 59.249 


285 


122 


258 


22.88 




Q 




J2347-1856 


23 47 08.6267 


-18 56 18.858 


244 


-205 


-134 


22.68 




G 





References. (1) McCarthy et al. (1996); (2) Healey et al. (2008); (3) Hall et al. (2002); (4) Gregg et al. (2002); (5) Ahn et al. (2012); (6) Tonry & 
Kochanek (1999); (7) Lawrence et al. (1995); (8) Wright & Otrupcek (1990); (9) Hook et al. (1996); (10) Vermeulen & Taylor (1995); (11) Falco 
et al. (1998); (12) Bade et al. (1998); (13) Carilli et al. (1998); (14) Hewitt & Burbidge (1993); (15) Fassnacht & Cohen (1998). 



Article number, page 10 of 10 



